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Abstract 

In random cellular systems, both observation and maximum entropy inference give 
a specific form to the topological pair correlation: it is bi-afhne in the cells num- 
ber of edges with coefficients depending on the distance between the two cells of 
the pair. Assuming this form for the pair correlations, we make explicit the condi- 
tions of statistical independence at large distance. When, on average, the defects 
do not contribute, the layer population and the enclosed topological charge both 
increase polynomially with distance. In dimension 2, the exponent of the leading 
terms depend on sum rules satisfied, or not, by the maximum entropy coefficients. 
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1 Introduction 



Foams are random but the precise probability distributions describing their 
structure are still largely unknown. Leaving aside the geometrical and ener- 
getic details, obviously dependent on the system considered, we focus on the 
topological characteristics, believed to be more universal. Only a few inves- 
tigations have been devoted to correlations beyond nearest neighbours. Pair 
correlations at arbitrary distance were analysed by Fortes-Pina Rivier et 
al pi3] . Szeto et a/ [3], and measured in [1]. 
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Maximum entropy is one of the few methods able to predict some aspects of 
the probabihty distributions. Maximum entropy (maxent) arguments yield a 
specific form for the pair correlations at arbitrary distance 



Here, we examine the asymptotic behaviour of pair correlations and how the 
independence of occurrences at large distances constrains the parameters in 
the maxent formulae. 

The present analysis is mostly devoted to bi-dimensional foams, subject to ex- 
tensive theoretical and experimental investigations 
and where the progresses are ahead of those on 3D foams. 



2 Foam statistics 



A foam T divides space into = \T\ polygonal cells. Here T is viewed as a 
set of cells and \T\ represents the number of elements in the set. 



2.1 One cell statistics 



The local observable is n, the number of sides of each cell ( "polygonality"), or 
the topological charge q = 6—n. The fraction of n-sided cells is p{n) = N{n)/N. 

The generic and physically stable vertex coordination (or degree) is z = 3. 
Euler's theorem implies: 



(g) = (Q-n) = ^(6 - n) p{n) = 6x/N 



where y characterises the embedding space: x = 2 for the sphere, for the 
toruo- 



In the large N limit, (g) 0. The second moment is fj,2 = (q^) = {{n — 6)^). 




Fig. 1. A pair of cells, {k,n) = (6,5), at distance j. 




Fig. 2. Stratification of a Poisson-Voronoi froth into layers. 
2.2 Two cells, correlations 

The topological distance j between cells is measured in nearest neighbour 
steps. The jth layer around a given cell o, lay(j|o), is the set of cells at 
distance j from a (Fig. [2]). It has population Kj{o) = |lay(j|o)|. The average 
over n-sided central cells is {Kj{n)) and the overall average is {{Kj)). 

^ As stated here, equation ([1]) only holds if the foam fills the 2D space, without 
boundary. Otherwise, the edges on the boundary must be counted only once in 
{n)N = 2E — i^boundaryj implying a additional boundary term in ([1]) [21j. E is the 
total number of edges, internal and on the boundary. In particular, this is so for a 
bounded cluster in the plane, where x = 1- lu any reasonable case, {q) = (6 — n) ^ 
as A'' ^ oo. 
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The joint distribution p^p{k,n) — the probabihty that a (fc, ?T,)-sided pair of 
cells occur at mutual distance j — and the corresponding marginal distribution 
Sj{n) = J2kpf\k,'n) — the probability that a cell is at distance j from an 

n-sided on^ — satisfy Sj{n) = p{n). 

The correlator Aj{k,n) and correlation function gj{k,n) are defined by 



■9j{k,n) Sj{k) Sj{n). (3) 



p^/\k,n) = A,{k,n)^^^^j^ (2) 



Both account for the statistical dependence of the concurrent occurrence of 
a k and a n-sided cell at distance j. They only differ in the way they are 
normalised: the correlation function is 1 whereas the correlator is the mean 
population {{Kj)) in independent situations. The correlator is related to the 
layer population by [Il|2] 

Y^p{k)Mk,n) = {K,{n)). (4) 
it 



A similar identity follows from counting the edges ('polygonality' Pj of layer 
j) DP: 

Y^kp{k)Aj{k,n) = {P,{n)). (5) 
it 



The polygonality of a set is the sum of the individual polygonalities n. Because 
charges are additive, and si mp ly related to n by g = 6 — n, it is more natural 
to consider the total chargqlJ of layer j, Q (lay {j\o)) = QKj{o) — Pj{o). Here 
we are getting closer to [2J. Then averaging over n-sided central cells o yields, 
according to (jl]), ([5]), 



(g(lay(j|n))) =5:(6 - k) p{k) A^{k,n). (6) 
it 



3 Maximum entropy 



Maximum entropy arguments (maxent) and the recursion relation to be de- 
scribed in sec. [5] give the following expressions, also observed in numerical 



This is also the probability that a n-sided cell is at distance j from an other one. 
^ The charge of a set of cells is the sum of the individual charges q = 6 — n. 



4 



simulations, [2] 



Aj{k,n) = aj{k - 6){n - 6) + aj{n + k - 12) + bj, (7) 
{Kj{n)) = aj{n-6) + bj, (8) 

aj,aj,bj are real parameters for j = 1,2,.... The contributions from de- 
fects have been neglected. In the infinite foam limit, the average of ([H]) gives 
{{K,)) = 6,. 



4 Asymptotic freedom 

In normal systems of statistical physics, spatially distant events become un- 
correlated. In foams, this was first measured by [1]. As j — > oo, 

pf\k,n) Sj{k) Sj{n) ^ gj{k,n) 1. (9) 
With 



'^^'^''^-{K,ik)){K,in)) 

= ^ + {^-it^)')ik-6){n-6) (11) 

and maxent, asymptotic de-correlation holds if and only if 



A sufficient condition is that the ratios ^ and P- ^ 0, meaning that, in 
both Aj and Kj, the n dependent terms would be dominated by the constant 
one (bj) at large j, a sensible result. In this limit, ([H]) implies {Kj{n)) — >■ {{Kj)) 
and then Sj{n) — > p{n)] the effect of the central condition vanishes, another 
manifestation of asymptotic de-correlation. In (Q, sj can then be replaced by 
p. 



5 Recursion equation 

The following equation was derived in [1] and in [2|22] : 

A{K,in))+ ^qM> {Kjin)) = ^ 0. (13) 
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The curly bracket quantities -< qjiji) >- (resp. -< rrijin) >~) are the cellu- 
lar charge (resp. sidedness) averaged over the neighbours at topological dis- 
tance i to an n-side cell. They are defined by -< qj{n) >- = 6— -< rrijin) >-= 
((5(lay(j|n)))/ {Kj{n)). AKj = Kj^i — 2Kj + Kj_i is the discrete laplacian. 



The right hand side is due to the presence of defects, the cells of layer j which 
have no edge in common with the next layer, j + 1 (Fig. [3]). It is assumed that, 
on average, this contribution vanishes: {Ij{n)) = 0. 




Fig. 3. Defects (left) and an inclusion (right), marked hatched. 



Using (ini) for the total charge (Q(lay(j|r;,))) =-<qj{n)>- {Kj{n)) allows to put 
the recursion relation ( JT3l) in the following form: 



A(ir,(n)) = Y.{k - 6) p{k) Aj{k, n). (14) 

k 



6 Asymptotic behaviour and sum rules 



With the maxent form ([7]), ([8]) of the correlator Aj{k, n) and population {Kj{n)), 
the recursion relation implies the following system, with aj = aj and jSj = 
bj - 6a f 



A 



a,- 











= , 









(15) 
(16) 
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The initial conditions follow from {Ko{n)) = and {Ki{n)) = n; fi2 was 
defined in sec. 12.11 

The solutions involve fk = — /U20"fc, conjectured to decrease fast at large k. 
Then, J2k=L fk — ^ 5* as j — oo, so that f/S defines a normalised distribu- 
tion (like probabilities except that some fk may be negative). Let M„ be the 
moments of this distribution: 

S = J2fk, M^ = Wkyk. (17) 

k>l ^ k>l 



The solutions to (fT5|) . (fT6!) are 



«, = E&/,+j(i-Ei=U), 

= ^ [(j + l)j(j - l)(l - Efk) + {3f - l)Ekfk 

- SjEk'fk + Ek'fk] + 6E(J - k)fk, 



(19) 



where the sums run from = 1 to j — 1. As j — > oo, they behave like 



A- 



■j(l-5) + 5Mi, 
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S 



(j + l)j(j-l)(l/5-l) + (3f -l)Mi 



-3jM2 + M3]+65(j-Mi; 



(20) 



(21) 



The leading order will be lower if the {fk = —fJ'2<^k} satisfy certain sum rules. 
Indeed, if 5 = 1, the estimates become 



A- 



Ml 
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{3f - l)Mi - 3jM2 + Mg 
-6(j-Mi), 



f22) 



(23) 



so that the average population {Kj{n)) would not grow faster than j^. 



The so-called Euclidean scaling [3f23] . {Kj{n)) ~ j in 2D, requires a second 
condition: 

Mi = -i^2Y.k(Tk = 0. 

k>l 
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7 Topological charge 



The purpose of this section is to relate the moments S and to the average 
charge contained in a ball or in a layer. It does not enforce the precise values 
found in the previous section, as we first hoped; but it shows that the sum 
rules and asymptotes, for charges and populations, are compatible. 



7.1 Charge enclosed in a ball 



Let Bj{o) be the union of the cells at distance at most j from an arbitrarily 
chosen cell o. The total charge enclosed in the cluster Bj{o) satisfies, averaging 
over n-sided origins o. 



(Q(S,)) = (6-n-/i2^ 



k=l 



+ 



/i2 



k=l k=l 



(24) 



where {Q{Bj)) stands for {Q{Bj{n))). This expression follows from inserting 
the solutions aj,Pj into ([71 [8]) and either summing the layer charges from 
to j or integrating the difference equation f|T3l) . It is consistent with fl2T]) . 
Indeed, if the individual charges are normal random variables, averaging to 
zero by Euler equation ([T]), their sum should behave like fluctuations: 



|(Q(S,(n)))|< const {\B,{n)\y/'. 



(25) 



Now, if 5 7^ 1, the population {\Bj{n)\) oc Eto(^K^)) ~ f and {Q{Bj{n))) 



J 



as in ( !24l) without further condition. 



If, on the other hand, the charge fluctuations are less, {Q{Bj{n))) = o{j 
then S = 1 must hold, predicting a layer population ~ j^, according to 



A possible interpretation of S" = 1 can be deduced from (HM . Define the 
average {{Qj)) = J2nP{^){Q{Bj{n))) . It coincides with the mean charge of a 
ball centred at a neutral cell: {Q{Bj{6))) . Then, using fl2^ . we can calculate 
the mean excess charge in the ball due to conditioning on n-sided central cells: 



{Q(B,{n)))-{{Q,)) = {6-n) 



E fk 

k=l 



(26) 
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As j oo, this excess charge becomes (6 — n)(l — S*); then 5* = 1 means global 
neutrality: the central charge q = 6— n is exactly screened by the excess charge 
(6 — 72)5' in the layers around. 



7.2 Layer charge and Aboav-Weaire's law 



Inserting the maxent formula ([8]) for {Kj{n)) into f|T3|) and using (|T5|) gives 
an expression for the average charge in layer j: 

(Q(lay(j|n))) = {n - 6)/, - fi2a„ (27) 



or, dividing by the mean layer population, 

^g,(n)^=M^^ = i!^^M:-/^ (28) 



This is a fractional linear function of n as in Aboav-Weaire's law [2lj, which 
is in fact the first case, j = 1, of this sequence labelled by j. Because double 
contacts are negligible in foams, the first layer population is just the number 
if sides of the central cell: Kxip) = {Ki{n)) = n and -< gi(n) >~= {qi{n)) = 
6 — (mi(n)), where (mi(n)) is the mean number of sides (the polygonality) 
of the first neighbours of n-sided cells. So, with ( lT6l) . equation ( l28i) for j = I 
specialises to 

{q,{n)) = ~ ' = /, - ^J^. (29) 
n n 



Compared to the usual form of Aboav-Weaire's law plp5] : (mi(n)) = 6 — 
a + (6a + ^2)/^, eq. (129|) gives an interpretation of the parameter a as the n k 
coefficient in the pair correlator: a = /i = —^2<^i- 



8 Conclusion 



To summarise, the major part of this article is a brief review of what is known, 
so far, on correlations in foams beyond first neighbours. This includes the def- 
initions and sum rules in sec 2, the (bi)affine form of the correlations in sec. 3, 
the recursion equation in sec. 5 and the solutions in sec. 6. Except for the affine 
ansatz, controversially [25j justified by maximum entropy arguments [2|5] . and 
the related hypothesis that defects do not contribute on average, which we as- 
sumed from the start, our purpose has been to restrain approximations or ad 
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hoc substitutions to a minimum. Some approximations (factorisations, etc.) 
give interesting perspectives and will be treated in an other article [26] . 

The new results are mainly contained in sections HI [6l [71 First, we have shown 
that the maxent affine ansatz is compatible with asymptotic freedom, and how 
the correlation decay affects the affine coefficients. 

Next, the decay of the correlations, or the rate of increase of the layer popula- 
tions, at large distance are tightly related to sum rules for the coefficients (Tj in 
equ. ([71), specifying the screening of (topological) charges. Charge neutrality 
means that a given charge is surrounded by a cloud of total opposite charge. 
In electrostatics, charge neutrality results from energy bounds and shielding 
the long range Coulomb field. This strong Debye screening implies subnormal 
charge fluctuations [27] . 

In foams, as far as we can see, no electric field constrains the topological 
charge. Global neutrality is a consequence of Euler equation Neutrality 
seems also true at an intermediate scale. The charge fluctuations result from 
a strange compromise between statistical disorder and geometrical constraints 
[2S] . In sec. [HI the order of the asymptotic polynomial behaviour of the layer 
population is related to specific sum rules satisfied by the correlation coeffi- 
cients. In turn, these sum rules command the overall charge (fluctuations) in 
large domains. 

The conclusion is that both populations and charges are consistently related. 
The next question is what do we need to get more precise, or more predictive, 
estimates. 
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